I 

I 


AD-A070  990  RENSSELAER  POLYTECHNIC  INST  TROY  N Y DEPT  OF  MECHANI— ETC  F/G  20/12 

THE  GROWTH  OF  WAVE  DISCONTINUITIES  IN  PIE2ELECTRIC  SEMICONDUCTO— ETC (U) 
JUN  79  M F MCCARTHY » H F TIERSTEN  N00014-76-C-0368 

UNCLASSIFIED  TR-28  NL 


DDC  FILE_COPY  M1A070990 


Rensselaer  Polytechnic  Institute 
Troyf  New  York  12181 


THE  GROWTH  OF  WAVE  DISCONTINUITIES  IN 
PIEZOELECTRIC  SOlICONDUCTORS 


by 


M.F.  McCarthy  and  H.F.  Tiersten 


Office  of  Naval  Research 
Contract  N00014-76-C-0368 
Project  NR  318-009 
Technical  Report  No.  28 


June  1979 


DDC 


7 

j|r?rrM?nn  nr? 

l 

f 

'<  JUL  11  1979 

J 

IktsEffTrEi 

yj 

D 


Distribution  of  this  document  is  unlimited.  Reproduction 
in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 

DISTRIBUTION  STATEMENT  'K 

Approved  for  public  releaso; 

Distribution  Unlimited 


i 


Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 


THE  GROWTH  OF  WAVE  DISCONTINUITIES  IN 
PIEZOELECTRIC  SEMICONDUCTORS 

i 

by 

M.F.  McCarthy  and  H.F.  Tlersten 


Accession  For 


NT IS  GKA4I 
DDC  TAB 
Unannounced 
Justification 


7 


1 button/ 


L. 


Mllty  Codes 

Avail  and/or 
special 


Office  of  Naval  Research 
Contract  N00014- 76-C-0368 
Project  NR  318-009 
Technical  Report  No.  28 


June  1979 


DDC 

[n^E£3EJlCl/7IEp 

JUL  11  I9T9 


Distribution  of  this  document  is  unlimited.  Reproduction 
in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 


jTSEQTTE 

D 


DISTRIBUTION  STATEMKNT  n 

Approved  for  public  releaso; 
Distribution  Unlimited 


1 


Unclassified 


SECURITY  CLASSIFICATION  OF  This  PAOE  fWA on  Data  Bn lorod) 


"4 


| REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER 

2.  OOVT  ACCESSION  NO. 

I.  RECIPIENT'S  CATALOG  NUMBER 

No.  28 

1.  TIHI  IKAUHN.) 7 -- 

THE  GROWTH  OF  ^AVE  DISCONTINUITIES  IN  / /2f 

) PIEZOELECTRIC  SEMICONDUCTORS  / ( ( > ^ 

.8.  TYPE  QF  REPORT  A BER.Op  COVERED 

) Technical  jtepmt  « ^ / 

~ — 

7.  AUTHORfal.. 

M.F. /McCarthy 


^McCarthy 


Tiers ten 


PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Department  of  Mechanical  Engineering, 
Aeronautical  Engineering  & Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 
Physics  Branch 


M MONITORING  AGENCY  NAME  S AOORESSflf  dllloront  tram  Controlling  Olll CO) 

(gjirg 


S CON' 


CONTRACT  OR  GRANT  NUMBERS «) 


N^0014-76-C-^368a 

"'N<.*-EN(rnio-M 

TtHPOIYRM  BLBMeHPl  PRSaBBWi  T AIK  -- 


AREA  • WORK  UNIT  NUMBERS 

NR  318-009 


V'j. 


32 


IS.  SECURITY  CLASS,  (ot  tfilo  roport) 

Unclassified 


ISa.  DECLASSIFICATION/ DOWNGRADING 
. SCHEDULE 


IS  DISTRIBUTION  STATEMENT  (ot  thlo  Roport) 

Distribution  of  this  document  is  unlimited. 


DISTRIBUTION  STATEMENT  fl. 

Approved  for  public  release; 
Distribution  Unlimited 


IT.  DISTRIBUTION  STATEMENT  (ot  tho  obotrmc t ontorod  In  Block  30.  It  dltloronl  troai  Roport) 


IS.  SUPPLEMENTARY  notes 


19.  KEY  WOROS  (Continue  on  rereree  aide  II  neceeeary  and  Identity  by  block  number) 


Electroelasticity 

Elasticity 

Nonlinear  Interactions 
Wave  Propagation 
Weak  Waves 

Induced  Discontinuities 


Wave  Growth 
Singular  Surfaces 


Deformable  Semiconductors 
Piezoelectric  Semiconductors 
Wave  Discontinuities 
Acoustoelectric  Domains 
Acceleration  Waves 
l^ySemiconduc tion 

*10.  (Continue  on  rereree  aldo  It  nocommmry  and  Identity  by  block  number) 

^ The  reference  coordinate  description  of  the  general  nonlinear  differ- 
ential equations  describing  the  interaction  of  finitely  deformable,  polar- 
izable, intrinsic  n-type  semiconductors  with  the  quasi-static  electric  field 
is  applied  in  the  study  of  acceleration  waves  in  piezoelectric  semiconductors 
As  a consequence,  the  mechanical  and  dielectric  nonlinearities  are  included 
in  the  treatment  as  well  as  the  semiconduction  nonlinearity.  The  general 
equation  for  the  propagation  velocity  of  the  disturbance  is  obtained  as  a 


DO 


FORM  im 
I JAN  71  Wl 


A> 


EDITION  OF  I NOV  BE  IS  OBSOLETE 
S/N  0 103*014*  SS01  I 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAOE  (Wbon  Dim 


fOI 


V- 


o 


Unc last! fled 


-H-UNITY  CLASSIFICATION  OF  THIS  PAlIi*Ii«t  Dm*  SkimmI) 


"^function  of  the  atate  of  the  material  Inmed lately  ahead  of  the  wavefront.  In 
Che  special  case  of  plane  waves  entering  a homogeneous  sCeady  state,  the 
growth  equation  for  the  amplitude  of  the  acceleration  wave  la  determined  and, 
of  course,  the  propagation  velocity  and  coefficients  In  the  growth  equation 
depend  on  the  propagation  direction,  but  otherwise  are  constant.  The  relation 
between  acceleration  waves  and  the  formation  and  propagation  of  acousto- 
electric  domains  is  indicated.  The  solutions  of  the  growth  equation  indicate 
the  formation  of  a shock  in  a finite  time  for  conditions  conducive  to  domain 
formation  except  in  certain  unusual  cases  possibly  occurring  with  purely 
transverse  acceleration  waves.  In  the  course  of  the  treatment  the  condition 
for  the  threshold  field  for  domain  formation  Is  determined  under  quite  general 
circumstances. m When  the  electrical  conduction  equation,  which  can  be  quite 
general  in  thlsT  treatment,  is  specialised  to  the  simple  form  usually  employed 
for  anisotropic\ semiconductors,  the  aforementioned  more  general  condition 
reduces  to  the  anisotropic  generalisation  of  the  well-known  elementary  result. 
In  addition,  the\ behavior  of  weak  waves  is  discussed. 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  FA0«P**WI  DM*  Sxt*r*« 


THE  GROWTH  OF  WAVE  DISCONTINUITIES  IN  PIEZOELECTRIC  SEMICONDUCTORS 


M.F.  McCarthy 

National  University  of  Ireland 
University  College 
Galway,  Ireland 

H.F.  Tiersten 

Department  of  Mechanical  Engineering, 
Aeronautical  Engineering  & Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 


ABSTRACT 

The  reference  coordinate  description  of  the  general  nonlinear  differ- 
ential equations  describing  the  interaction  of  finitely  deformable,  polar- 
izable, intrinsic  n-type  semiconductors  with  the  quasi-static  electric 
field  is  applied  in  the  study  of  acceleration  waves  in  piezoelectric  semi- 
conductors. As  a consequence,  the  mechanical  and  dielectric  nonlinearities 
are  included  in  the  treatment  as  well  as  the  semiconduction  nonlinearity. 

The  general  equation  for  the  propagation  velocity  of  the  disturbance  is 
obtained  as  a function  of  the  state  of  the  material  immediately  ahead  of 
the  wavefront.  In  the  special  case  of  plane  waves  entering  a homogeneous 
steady  state,  the  growth  equation  for  the  amplitude  of  the  acceleration 
wave  is  determined  and,  of  course,  the  propagation  velocity  and  coefficients 
in  the  growth  equation  depend  on  the  propagation  direction,  but  otherwise  are 
constant.  The  relation  between  acceleration  waves  and  the  formation  and 
propagation  of  acoustoelectric  domains  is  indicated.  The  solutions  of  the 
growth  equation  indicate  the  formation  of  a shock  in  a finite  time  for 
conditions  conducive  to  domain  formation  except  in  certain  unusual  cases 
possibly  occurring  with  purely  transverse  acceleration  waves.  In  the 
course  of  the  treatment  the  condition  for  the  threshold  field  for  domain 
formation  is  determined  under  quite  general  circumstances.  When  the  elec- 
trical conduction  equation,  which  can  be  quite  general  in  this  treatment, 
is  specialized  to  the  simple  form  usually  employed  for  anisotropic  semi- 
conductors, the  aforementioned  more  general  condition  reduces  to  the  aniso- 
tropic generalization  of  the  well-known  elementary  result.  In  addition, 
the  behavior  of  weak  waves  is  discussed. 
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1.  Introduction 

In  a previous  investigation  the  theory  of  one-dimensional  acceleration 
waves  was  applied1  to  a one-dimensional  version  of  general  rotationally  in- 
variant nonlinear  electroelastic  equations  derived  earlier  from  a well- 

2 1 
defined  macroscopic  model  of  deformable  semiconductors.  In  that  treatment 

an  analytical  description  of  the  formation  and  propagation  of  purely  longi- 
tudinal acoustoelectric  domains  in  piezoelectric  semiconductors  was  obtained. 
The  analysis  indicated  that  for  electric  fields  above  a threshold  value  the 

amplitude  of  the  acceleration  wave  would  always  increase  without  bound  and  become 
a shock.  A natural  and  logical  extension  of  the  previous  one-dimensional 
work  is  the  treatment  of  three-dimensional  acceleration  waves,  in  which 
acoustoelectric  domains  with  transverse  mechanical  displacement  components 
can  be  considered.  Recently,  in  the  case  of  the  quasi-static  electric  field 

the  general  nonlinear  electroelastic  equations  for  deformable  intrinsic 

2 3 

n-type  semiconductors  were  transformed  from  the  unknown  present  coordinate 
description  to  the  known  reference  coordinate  description,  which  is  the  form 
needed  here  and  in  general  for  the  treatment  of  problems. 

4-12 

In  this  paper  the  theory  of  three-dimensional  acceleration  waves 

3 

is  applied  to  the  above-mentioned  reference  coordinate  description  of  the 
general  rotationally  invariant  nonlinear  electroelastic  equations  for  de- 
formable intrinsic  n-type  semiconductors  in  order  to  analytically  describe 
the  formation  and  propagation  of  acoustoelectric  domains,  with  both  trans- 
verse and  longitudinal  components  of  mechanical  displacement,  in  piezoelectric 
semiconductors  subject  to  high  electric  fields.  The  analysis  results  in  an 
expression  for  the  amplitude  of  the  acceleration  wave  (or  domain)  which 
exhibits  the  competition  between  dissipation  due  to  electrical  conduction 


i 


i 


2. 


and  the  semiconduction  and  mechanical  nonlinearities  in  producing  decay  or 
growth  of  the  acceleration  wave  (or  domain).  As  in  the  case  of  the  purely 
longitudinal  acceleration  wave  treated  earlier1,  the  possibility  of  the 
amplitude  of  the  more  general  three-dimensional  plane  acceleration  wave 
increasing  without  bound  and  becoming  a shock  is  clearly  indicated.  How- 
ever, in  the  special  case  of  purely  transverse  acceleration  waves,  circum- 
stances can  exist  under  which  it  is  not  possible  for  the  amplitude  to  grow. 
Nevertheless,  if  any  longitudinal  motion  is  present  in  the  acceleration 
wave,  the  possibility  of  the  amplitude  increasing  without  bound  always 
exists.  During  the  course  of  the  analysis  the  expression  for  the  velocity 
of  the  wave  (or  domain)  as  a function  of  the  state  of  the  material  immedi- 
ately ahead  of  the  wavefront  naturally  is  obtained. 

2.  Fasic  Formulae  and  Equations 

The  macroscopic  model  of  an  elastic  intrinsic  n-type  semiconductor 
employed  in  Ref . 2 consists  of  three  interacting,  interpenetrating  continue, 
which  consist  of  (i)  a lattice  continuum  which  has  a positive  charge  densityj 
(ii)  a bound  electronic  continuum  which  has  a negative  charge  density  ami 
which  can  displace  slightly  from  the  lattice  continuum  and  thus  produce 
electric  polarization,  and  (iii)  a free  electronic  continuum  which  has  a 
negative  charge  density,  negligible  inertia  and  is  a conducting  compressible 
fluid  which  experiences  a force  of  resistance  from  its  motion  with  respect 
to  the  lattice  continuum. 

Initially,  the  lattice  continuum  and  the  bound  electronic  continuum 
all  occupy  the  same  region  of  space  and,  hence,  have  the  same  reference 
coordinates  X^.  The  motion  of  a point  of  the  lattice  continuum  is  described 
by  the  mapping 


3. 


yi-yi(Vt) > (2.1) 

which  is  one-to-one  and  differentiable  as  often  as  required.  Here  the  y^^ 
denote  the  present  coordinates  of  material  (lattice  continuum)  points  and 
X^,  the  reference  coordinates,  and  t denotes  the  time.  We  consistently  use 
the  convention  that  capital  indices  denote  the  Cartesian  components  of  X 
and  lower  case  inices,  the  Cartesian  components  of  £.  A comma  followed  by 
an  index  denotes  partial  differentiation  with  respect  to  a coordinate 


9 

> 


i 


(XR,t) 


f 


(2.2) 


and  the  summation  convention  for  repeated  tensor  indices  is  employed. 

Since  reference  coordinates  are  employed  in  our  study  of  the  propaga- 
tion of  acceleration  waves  in  elastic  semiconductors,  the  integral  forms 
of  the  equations  required  in  this  work  consist  of  Eqs. (2.41)  - (2.44)  of 
Ref. 3,  Which  we  reproduce  here  in  the  form 


I ',L<KLj+MLj-6?j,aSo-£  J "oW 

S J j j v 

o o 

J ffaV 

S V 

o o 

J WLdVo“J  NL(CP  + Cf>e)  dSo  , 


(2.3) 

(2.4) 

(2.5) 


J "A".- 

s 


(2.6) 


o o 

where  N denotes  the  outwardly  directed  unit  normal  to  a reference  element 
If 

of  area  and  S denotes  the  surface  enclosing  the  reference  volume  V . 

o o 

Equations  (2.3)-  (2.6)  are  the  reference  integral  forms  of  the  conservation 
of  linear  momentum  of  the  combined  continuum,  the  charge  equation  of 


i 

; 


inn^rtfwi 
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electrostatics,  the  conservation  of  linear  momentum  of  the  free  electronic 
contimnan  and  the  conservation  of  total  electric  charge.  In  Eqs.  (2.3)  - (2.6) 
V denote  the  reduced  mechanical  Eiola-Kirchhof f streps  tensor, 


the  reference  free-space  Maxwell  electrostatic  stress  tensor  and  the  refer- 

e 


ence  free  electronic  pressure  tensor,  respectively;  v , 3 , w and  Q denote 

j L L L 

the  velocity  of  the  solid,  the  reference  electric  displacement  vector,  the 


reference  local  electric  field  exerted  on  the  free  electronic  fluid  and  the 


reference  electric  current  vector,  respectively;  pQ,  n,  cp  and  cp  denote  the 


reference  mass  density,  net  reference  charge  density,  electric  potential  and 
free  electronic  chemical  potential,  respectively;  and  d/dt  is  the  material 
time  derivative.  The  associated  constitutive  equations  and  additional  re- 
quired relations  take  the  form 


» MLj"JXL,iTij 
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SX 


* . ■ JX  p 3 m e JX  . E . - p , 

Lj  L,  j * L o L,  l i ro  do  1 


e^e 


p'-<pe>2^, 

an 


(2.7) 


ES  1 1 i 


4i  " ” i ’ EKL  " 2 (CKL  ” ’ CKL*yi,K?i,L* 


(2.8) 


x"x(lW,i)  » «e-«e»e).  9e-aft*e«e)/a»e, 


(2.9) 


_ES 


where  T. E and  C denote  the  free-space  Maxwell  electrostatic  stress 
ij  KXj  KXj 

tensor,  the  material  (reference  or  Lagrangian)  strain  tensor  and  Green' s 

0 

deformation  tensor,  respectively;  E^,  <$L  and  denote  the  Maxwell  electric 
field,  the  reference  (or  rotationally  invariant)  measure  of  the  electric 
field  and  the  rotationally  invariant  constitutive  vector  that  accounts  for 


li 
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the  relative  flow  velocity  of  the  free  electronic  fluid,  respectively;  p , 


X,  p®  and  «*  denote  the  free  electronic  pressure,  a particularly  convenient 


thermodynamic  state  function  related  to  the  stored  internal  energy  per  unit 
mass  of  the  deformable  solid,  the  free  electronic  charge  density  and  the 
stored  internal  energy  per  unit  charge  of  the  free  electronic  fluid,  respect- 


ively and  c is  the  permittivity  of  free  space, 
o 


When  the  variables  are  appropriately  differentiable,  from  (2.3)-  (2.6), 
we  obtain  the  differential  equations 


where 


S ■ Q v 

Lj,L  PoV 

(2.10) 

(2.11) 

(2.12) 

A. 

*1,1*“-° 

(2.13) 

(2.14) 

— e r 
P-JP  +oP  , 


and  we  have  employed  the  dot  notation  for  partial  differentiation  with 
respect  to  time.  We  now  note  that  we  have  an  additional  relation  between 
the  net  reference  charge  density  p and  the  free  electronic  charge  density  pe, 
which  can  be  written  in  the  form 

(2.15) 

where  qP~  is  the  reference  residual  lattice  charge  density,  which  is  a 
constant.  Prom  (2.7)^.,  (2 . S) 7,  <2-9>2-4  an<J  we  can  v*ite 

(2‘16) 

where 

V"**  • 
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defi«e  x by13 
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thM  * ««„e  „ (2  ej 
— — r •9)— 
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We  can  write  ^ " \i\i> 


„ X*X(E  , 

(2.7,  „ 


f2.29; 


(2.20) 


dj  a 

a ,• , -Jki 

and  the  chain  rule  of  diff  <yi,M)  ' 
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7t  r 


(2.21) 
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f . i df/dt 

i ” If  ,.l  ’ Un""  If  . I * 


Corresponding  to  a,  we  have  the  alternative  representation  of  the  surface  £ 
in  (X, t)  space  by  means  of  the  equation 

3(X,t)  - f (y(X,t),t)  - 0 . (3.3) 

The  unit  vector  N normal  to  £ and  its  speed  of  propagation  are  given  by 


Nk  = 


> un“- 


It  is  a simple  matter  using  (3.2)-  (3.4)  and  the  chain  rule  of  different- 


iation to  show  that 


F..  n. 
ik  l 


U„=»U 


K ’ N ip.Kl 


“'I'V’iJ  - 


(3.5) 


Where 


PiK  yi,K  > 


U = u - v • n . 

n ~ 'w  * 


is  the  local  speed  of  propagation  of  the  surface. 

Let  |(X, t)  be  a function  Which  suffers  a jump  discontinuity  across 
the  surface  £,  but  is  a continuous  function  everywhere  else.  We  define 
the  jump  [+]  in  the  function  A to  be 

m = f-*+,  <2 


where  t and  f are  the  limiting  values  of  t immediately  behind  and  just 

in  front  of  a point  lying  on  £,  respectively.  The  surface  £ is  said  to 

be  an  acceleration  wave  if  the  fields  y. (X, t),  v.  (X.t)  and  F.,(X, t)  are 

1 ~ 7 l lL 

continuous  everywhere  but  <r.  (X.t).  F.(X,t)  and  F..  ..(X.t),  as  well  as  all 

1 ~ 7 1I<  ~ 1L,  M ~ 7 


8. 


higher  order  partial  derivatives  of  y^(X, t),  suffer  jump  discontinuities 

across  E,  but  are  continuous  functions  everywhere  else.  From  the  geometric 

14 

conditions  of  compatibility  on  the  jump  in  the  gradient  of  a continuous 
function  and  the  kinematic  condition  of  compatibility15,  we  can  obtain1** 


r.  k in, 


JPiK,li-WL-aiFpKVVV  V 

VA-  °*iW 

• 2 2 2 2 
[v.)«U„s  = U a , , s . = B a , B =F.  F,  n.n,  . 
~ i~  N i l ' l n i’  n lK  jK  l j 


(3.8) 


At  this  point  it  should  be  noted  that  we  need  make  no  assumptions  with  regard 

to  the  continuity  properties  enjoyed  by  the  electric  potential  cp(X, t)  = 

0 

cp(y(X,t),t)  or  the  free  electronic  charge  density  p.  (X,t)  apart  from  assuming 
that  at  points  not  on  £ these  functions  together  with  their  partial  deriva- 
tives of  all  orders  are  continuous.  The  vector  a is  called  the  amplitude 
vector  of  the  acceleration  wave.  If  we  write  a = ar,  where  r'niO.  I r|  = 1, 
then  if  a>0  the  wave  is  said  to  be  expansive.  While  a wave  for  which  a<0 
is  said  to  be  compressive.  If  r ** n the  wave  is  longitudinal,  while  if  r • n * 0 
it  is  transverse. 

The  jump  conditions  across  a surface  of  discontinuity  E can  readily  be 
obtained  from  the  integral  forms  in  (2.3)-  (2.6)  along  with  the  fact  that 
remains  bounded.  The  resulting  jump  conditions  thus  obtained  consist  of 
Eqs.  (2.46)1_2,  (2.47)  , (2.49)  and  (2.50)  of  Ref. 3,  which  are  required  in 
this  work  and  we  reproduce  here  in  the  form 


V4cL  -Vffi-0’ 

(3.9) 

W * Wvj!-°- 

(3.10) 

iA-°.  m-°, 

(3.11) 

■UM  - 0 • 

(3.12) 
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In  view  of  (2.9)^,  (3.11)^  may  be  written 
• e e-  «e  e+ 

9 (p,  ) - cp  fti  ) - 0 , (3.13) 

and  if  we  assume  that  dcpe  (pie)/c!p,e  / 0 it  follows  (the  argument  is  given  by 
6 0 

Coleman  and  Gurtin  ) that  p is  continuous  across  E,  i.e., 

D**]  - 0.  (3.14) 

Furthermore,  since  (2.15)  is  of  the  form 

S - Jp.®  + JJ.r  , o^T  - constant  (3.15) 

it  follows,  since  J=det  F^K  is  continuous  across  E,  that 

rtE]  -0.  (3.16) 

Next,  since  cp  is  continuous  across  E,from  the  geometric  condition  of  com- 
14 

patibility  and  (2 . 8 ) 7,  we  have 

Thus,  in  view  of  (2.23)2,  we  may  write  (3.12)  in  the  form 

NA<ERL’<ti!V  -°.  <3-ie> 

and  if  we  assume  that  d^/55^/0  it  follows  from  (3.18)  that  5=0  and,  hence, 
from  (3.17)  that 

I^kL"0’  (3.19) 

and  thus  the  reference  electric  field  S is  continuous  across  E.  Now.  in 

»v  ' 

view  of  (3.14),  from  (2.17)  and  the  geometric  condition  of  compatibility 
we  have 

Igk l-%,  w* iX>  (3'20) 

so  that,  on  using  (2.16),  (3.14),  (3.16),  (3.19)  and  (3.20),  we  may  rewrite 
(3.9)  in  the  form 


10. 


Prom  the  assumption  that  d^^/SG^O  follows  from  (3.21)  that  <u  - 0 
so  that 

JG^-O,  (3.22) 

which  means  that  is  continuous  across  E.  As  a consequence  of  (2.16), 
(2.8)g_^,  (2.23),  (3.14),  (3.19),  (3.22)  and  the  definition  of  an  accelera- 


tion wave,  we  have 


JV’Ai-Mci-0- 


(3.23) 


and  thus  Eqs.  (3.9)-  (3.12)  are  satisfied  identically.  Other  important  con- 
ditions resulting  from  (3.19),  (3.22)  and  the  conditions  of  compatibility14-1*’ 

are 

I*K,  li  ' " ’ I*  kJ.' - - °Vl/  I\i  * V*K  ’ “ ‘ I W,  RS~  ’ 

W "-iW^nsl  * (3’24) 

Prom  (2.10),  (2.22),  (2.7)3  g,  (2.8)5_6,  (2.9)12  and  (2.21),  we  obtain 

(3.25) 


AjKpLFpL,K  + ^KlA^K  + HjKGK*Po'rj» 


where 


'jKpL  - SPpL  - -o  SFjKaFpL  - 3Pp[j  <JXK,))P  “‘C)  ' ‘jpS«L  + 


+ VpSC»LS‘ J«K,JXI.,p-XL,)XK,p>p0(pC»  > 


jKL 


L 

!> 

jK 


5e! 

an* 


and 


d2* 


A 2 A 

S m n _ _ m Q ^ X n m ^ '• 

KL  Mo  dE„T  • KNLS  H0  dE~3E77  ’ KNI.  Mo  501“  , 
KL  KN  Lo  L KN 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


11. 


Since  S is  by  hypothesis  a C -function,  it  follows  that  the  coeffi- 

cients  in  (3.26)-  (3.28)  are  continuous  across  the  singular  surface  Z. 

Thus,  taking  the  jump  in  Eq. (3.25)  across  Z and  employing  (3.8)^,  (3.22) 
and  (3.24)^  we  obtain 


VplWp  - “jklW  ».V)  • 


(3.30) 


We  now  need  the  expression  for  a in  terms  of  s.  Prom  (2.11)  and  (3.27)  we 


BiLKFiL,  K + *)aA,  K " ^ * 


where 


kkl  p< 

Li 


a2x 


(3.31) 


(3.32) 


Since  & (.,.)  is  a C -function,  it  follows  that  the  coefficients  in  (3.31) 
are  continuous  across  the  singular  surface  Z.  Thus,  on  taking  the  jump  in 
(3.31)  across  Z and  using  (3.8)^,  (3.16)  and  (3.24)^,  we  find 


where 


Cl.s  , 


Li  “ BiLKNLNK-  FiMBMLKNLNK  ’ <> 


(3.33) 


(3.34) 


The  substitution  of  (3.33)  into  (3.30)  with  the  aid  of  (3.26)  and  (3.34)^ 
yields  the  following  propagation  condition 


(3.35) 


where 


°jP- VP«A-V  -WkV^V  V>A"i  * 


FjLFpMCKLRMNKNR  + ^LjLp  ®jp 


(3.36) 


is  the  acoustic  tensor.  We  note  that  Q.  is  symmetric  and  for  fixed  n it  is 

jp 

a function  of  the  deformation  gradient  P , , the  reference  electric  field 

rir  k 

0 

at  the  wavefront,  but  is  independent  of  p , the  density  of  free  electronic 


charge  at  the  wavefront.  It  follows  from  (3.35)  that  the  amplitude  a of 


12. 


an  acceleration  wave  traveling  in  the  direction  n in  a piezoelectric  semi- 

A 

conductor  must  be  a proper  eigenvector  of  the  symmetric  acoustic  tensor 

2 

and  the  speed  of  propagation  UN  must  be  such  that  PQUN  is  the  corresponding 

A 

eigenvalue  of  Q . 

The  equation  (3.35)  was  derived  by  Truesdell17  for  acceleration  waves 

in  elastic  media.  It  has  since  been  derived  by  a number  of  authors  for 

7-9 

acceleration  waves  in  a variety  of  media  . We  note  in  particular  that 

the  acoustic  tensor  (3.36)  has  precisely  the  same  form  as  the  corresponding 

acoustic  tensor  which  occurs  in  the  theory  of  wave  propagation  in  elastic 
10 

dielectrics 

Equations  (3.35)  admit  a nontrivial  solution  if  and  only  if 


det  (Qjp-P0VjP)"°  * (3-37) 

and  this  equation  determines  the  possible  speeds  of  propagation  for  a given 
direction  of  propagation  n.  On  the  other  hand,  if  the  amplitude  a of  a wave 
is  known,  then  the  corresponding  speed  of  propagation  is  determined  by  the 
formula 

PoUn  " V*jap/aiai  • (3*38) 


Since  the  acoustic  tensor  $ is  symmetric,  it  has  three  real  eigenvalues. 
However,  at  this  stage  it  is  possible  that  all  of  these  eigenvalues  may  be 
negative  in  which  case  no  real  waves  will  exist  at  all.  We  now  wish  to 
record  the  conditions  which  guarantee  the  possible  existence  of  at  least 
some  real  waves.  A detailed  analysis  of  the  situation  for  purely  elastic 
materials  has  been  given  by  Truesdell18,  Truesdell  and  Noll19,  Wanq  and 
Truesdell20,  Chadwick  and  Currie21. 

Oice  the  deformation  gradient  and  electric  field  ahead  of  the  wave  are 
known,  it  follows  from  (3.5)^  (3.29),  (3.34),  (3.36),  (2.8)5_fi  and  (2.9^ 


that  the  acoustic  tensor  depends  on  n only.  That  is,  we  have 


VFriA-V  ■ VS’  - 


for  fixed  F and  <S  . If 
rL  K 


Qij (S>ninj  > 0 


(3.39) 


(3.40) 


for  all  unit  vectors  n,  the  material  may  be  said  to  have  positive  longitudinal 

•v7 
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piezoelectricity.  Truesdell  has  shown  that  when  (3.40)  is  satisfied  there 
exists  at  least  one  direction  in  which  a longitudinal  wave  may  exist  and 
propagate.  If  the  acoustic  tensor  is  strongly  elliptic  in  the  sense  that 


Qij j > 0 ’ (3-41) 

18 

for  all  unit  vectors  n and  [i.  it  then  follows  (Truesdell  ) that  there  is 
at  least  one  direction  of  propagation  in  which  a longitudinal  wave  and  two 
transverse  waves  with  orthogonal  amplitudes  may  exist  and  propagate.  In 
particular,  it  should  be  noted  that  if  the  strong  ellipticity  condition  (3.41) 
is  satisfied  and  if  the  deformation  and  electric  field  are  uniform  ahead  of 
the  wave  then  there  exists  at  least  one  direction  in  which  a plane  longi- 
tudinal and  two  plane  transverse  waves  may  exist  and  propagate  for  all  times. 
For  propagation  in  all  other  directions  under  the  above-mentioned  circum- 
stances the  three  plane  waves  are,  of  course,  not  necessarily  either  purely 
longitudinal  or  purely  transverse,  but  may  consist  of  an  admixture  of  all 
mechanical  displacement  components.  Nevertheless,  in  the  most  general  case 
if  the  propagation  velocities  are  distinct,  the  three  plane  waves  have 
mutually  orthogonal  mechanical  displacement  fields.  In  the  next  section, 
we  examine  the  manner  in  which  the  amplitudes  of  such  waves  vary  as  they 


traverse  the  material. 


14. 


4.  Growth  and  Decay  of  Plane  Acceleration  Waves 

In  this  section  we  derive  the  differential  equation  which  determines 
the  manner  in  which  the  amplitude  of  a plane  acceleration  wave  varies  as  it 
traverses  the  material.  It  is  assumed  that  the  material  ahead  of  the  plane 
wavefront  is  at  rest  in  a state  of  homogeneous  strain,  is  subject  to  a 
uniform  electric  field  and  that  the  charge  density  of  the  free  electronic 
fluid  is  uniform  and  constant  prior  to  the  arrival  of  the  wavefront. 

Hie  differentiation  of  the  equation  of  motion  (3.25)  with  respect  to  t, 
with  X fixed,  yields 


where 


with 


AjKpLFpL,  K + BjKL*L,  K + Hjx\  + ° j “ PoVj  * 


°j  " ^jKpI/jM^pL, K*qM  + BjKpLM^pL, K M 
+ H j KpLFpL, + B j KpML^L, KFpM 
+ BjKLM*L, K^M  + H jKpLGKFpL 


(4.1) 


+ W 


• e 


(4.2) 


2* 

is. 


IV  ■ . - SPc  +6Pc 

jKpLqM  3F  3P  jp  qN  KLMN  jq  pN  KMLN 

pL  qM 


+ 6pqF  jN°KNLM  + F jNFpTFqR°  KNLTMR  " ^M,  q (\,  jXL,  p 

- X _,X  )pe  + J (X  X .X  +X  .X  X 
L,  j K,p  F K,q  M,  j L,  p K,  j M,  p L,  q 


‘ XL,qXM,jXK,p-XL,jXK,qXM,p,P  ’ 


32s 


5J-  - Lb + F...P  .b. 


jKpLM  3Fr5J„  "jp“KLM  ' * jN‘pS~KNLSM’ 
pL  M 


.*6 


jKpL  *«3F 


- J(X  x - X X ) zz-  , 
K,j  L,p  L,j  K,p  ^e 


32s 


pL 


Kj  ~ 3 p 

BjKLM  " 3tfL3<fM  " FjNBNKLM  ’ Hjx““  JXK,  j ’ 


2*e 

3<pV 


i 


(4.3) 
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and 


'KLMNRT 


KNLSM 


bnklm 


8EM»dV 

Po5VS1^^ 

d2s 

KM 

r 

° dEKMdELSd<,M  ' 

d2S 

NK 

a3x 

4 W4»  ' 

(4.4) 


In  view  of  our  earlier  assumptions  on  the  continuity  of  the  response  func- 
tions the  coefficients  in  (4.3)  are  continuous  across  the  singular 

surface  E.  On  taking  the  jump  in  (4.1)  across  Z,  keeping  in  mind  our  assump- 
tions concerning  the  uniformity  of  the  rest  state  of  the  material  ahead  of 
the  wavefront  and  making  use  of  the  relation 

Ml.  = , (4.5) 


and  employing  (3.22)  we  obtain 

AjKpliFpL, kI  + BjKli^L, kI  + H jKi^xi  + 1C jl = PoIVjl 


(4.6) 


Where 


-Cj~  * AjKpLqMiFpL,  kI  + BjKpLM^IFpL,  kI  I^mJ> 


(4.7) 


Since  we  are  dealing  with  a plane  wavefront  entering  a uniform  state,  from 

the  compatibility  conditions  and  the  definition  of  an  acceleration  wave, 

^ 22 
we  have 


■FrL.Kl-WiC-  br-IWrM,sl  > 
+ TT- 

I^kI-^lV  s’-  IVs’rsI- 


(4.8) 
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Where  6f/6t  is  the  displacement  derivative  of  f.  Substituting  from  (3.8), 
(3.24)  and  (4.8)  into  (4.6),  we  obtain 

6s. 


"WV*L-  «V6p*(,iKLVL-  Wk-  si 


VjKPiVVVkV,  * 2D»<ffljKpLMVLnMSp  - V WWh-°-  ,4-9> 


Me  now  need  the  expression  for  cu  in  terms  of  s.  Tt>  this  end.  with  the 

'V  * 

aid  of  (2.16)  and  (2.8)5_^  we  rewrite  Eg. (2.13)  in  the  form 


^Krl/rL,  K + ^KL^L,  K + fW%  K + W + ^ ” ° 


(4.10) 


where 


. it , 

“ SFrL  * ^N^KLN'  EKMN  “ * EJCL  " &T  » 


^KL'ac"' 


(4.11) 


In  view  of  the  continuity  assumptions  on  ),  it  follows  that  the 

coefficients  in  (4.11)  are  continuous  across  the  singular  surface.  Hence, 

on  taking  the  jump  in  (4.10)  across  £ and  employing  (3.8),  (3.22)  and 

(3.24)  we  have 

l,4’ 

” ' v ‘-WWr  + “Wl  * I51  > 

where 


(4.12) 


V-  "VlW 


-1 


(4.13) 


Taking  the  material  time  derivative  of  (2.15),  we  obtain 


— • e • e 

p = J)i  +JP  , 


(4.14) 


the  jump  in  which,  with  <2.8)p  (3. 8)3,  (3.14),  (3.22),  the  kinematic  condi- 


15 

tion  of  compatibility  and  the  definition  of  an  acceleration  wave,  yields 


' 

J 
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where  we  have  employed  the  well-known  relation 

(4-16’ 

in  obtaining  (4.15).  The  substitution  of  (3.33)  and  (4.15)  into  (4.12) 
enables  us  to  write 


in  « vE  s , 
r r ’ 


where 


E - -E  nn-CE  nnl+ju  UCN  X 
r KrL  L K ^KLKLr  N K K,r 


(4.17) 


(4.18) 


We  now  must  express  3 in  terms  of  s and  b.  To  this  end  we  take 

<v 

the  material  time  derivative  of  (3.31)  to  obtain 


“BiLKFiL,  K * HKl\,  K " BiLrSKFrSFiL,  K 


-B.  (<$  F,  + F . S ) 

iLKM  M iL,K  iL  M,  K 


+ hklm<Vl,  K “ ^ 


where 


2 A 

\lm  " sjtst  " ~ po  ’ 

L M K L M 


(4.19) 


(4.20) 


which  is  continuous  across  the  singular  surface,  clearly,  all  the  coeffi- 
cients occurring  in  Eq.  (4. 19)  are  continuous  across  the  singular  surface  E. 
Thus,  on  taking  the  jump  in  (4.19)  across  E,  recalling  the  time  independent 
uniform  state  ahead  of  the  wavefront  and  employing  (4.15)  we  have 


SILkI^1L,k1  + kI  ~ PiLrSKi^rsIIPiL,  id 


* WV^L,a--Sl 


(4. 21) 


1 


Substituting  from  (3.8),  (3.24),  (3.33),  (3.34),  (4.8)  and  (4.15)  into  (4.2l), 


we  obtain 


where 


» - ' c vX,  iVi * Vi3  v j 


Mij  - <-Biio»c  * ♦ ('UiV™ 


(4.22) 


(4.23) 


We  now  substitute  from  (3.33),  (4.17)  and  (4.22)  into  (4.9)  and  employ 


(3.34)1  and  (3.36)  to  obtain 


where 


A 2.  2 j 

(Q,  -P  tJ„6 . )b  = 2p  u —r~  + a.  s +a.  s s , 
jp  o N 3p  p Ko  N 6t  jp  p jpq  p q ’ 


(4.24) 


(4.25) 


“jpq  UN^LjMpq  + UNfAjKpLqM  + 2^BjKpLMI'q+ ^ BjKLMLpLqlNKNLNM  ' t4-26) 


Note  that  Eqs. (4.24)  serve  to  determine  the  components  of  the  amplitude  b 
of  the  third  order  discontinuity  induced  by  the  acceleration  wave.  However, 
for  the  moment,  our  primary  objective  is  to  use  Eqs. (4.24)  to  obtain  the 
differential  equation  governing  the  evolutionary  behavior  of  the  amplitude  a 

<v 

of  the  acceleration  wave.  In  Sec. 6 we  study  the  solution  of  (4.24)  in 
greater  detail  and  we  discuss  both  induced  discontinuities  and  higher  order 


waves. 


Using  Eqs.  (3.5),  (3.8)4_g,  (4.23),  (4.25)  and  (4.26),  we  may  rewrite 


Eq.  (4.24)  in  the  form 


where 


a 2 2 -i 

(Q.  - p U 6.  )b  = 2p  0 — + o?.  a +or.  a a , 
jp  o N dp  p 'o  6t  jp  p DPq  P q' 


5jp*CU0HeLjnp  + vUHjKPsKnsZp, 


(4.27) 


(4.28) 


■ 


IL 
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5)Pl'"'I-JiiPq*U,,ljK«.WM*2BjKI,LML,*t  Wf  W.tVrVt ' <4”> 


and 


M » [-B 


pq  ®pKqLM  + 2^BpKLMLq  + ^ 'ScLM1  FrKFsLFtMnrnsnf 

We  may  write 

vav  y-i» 


(4.30) 


(4.31) 


A 2 

where  r is  the  unit  eigenvector  of  Q corresponding  to  the  eigenvalue  p UM, 


o N’ 


and  recall  that,  since  we  are  dealing  with  a plane  wavefront  propagating 
into  a uniform  region,  the  components  of  r are  constant.  Equation  (4.27) 


may  now  be  written 


l6Jp-  Po^e)p,bp-2|!o“2rj  I?  » <5)prp»»*  'Sjpq'pV*2 


(4.32) 


If  we  now  contract  (4.32)  with  r^  and  employ  (3.35)  and  (3.36),  we  find  that 
the  amplitude  a of  the  acceleration  wave  satisfies  the  equation 


where 


6a  0 2 

■jr—  »-u)a+p  a , 
6t  o o ’ 


a.  r r 

oi  - JP  3 P 

o 2 * 

2p  U 
o 

-5.  r.r  r 

b . - jpg 

° o r,2 

2p  0 

o 


(4.33) 


(4.34) 


(4.35) 


The  implications  of  Eq.  (4.3  3)  are  examined  in  detail  in  the  next  section. 
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5.  Implications  of  the  Growth  Equation 

Equation  (4.34),  which  governs  the  evolutionary  behavior  of  plane 
acceleration  waves  oriented  arbitrarily  and  propagating  into  a region  subject 


to  a time  independent  state  of  homogeneous  strain  and  uniform  electric  field 


in  a deformable  semiconductor,  is  of  the  same  form  as  Eq.  (5.6)  of  Ref.l. 

Clearly,  we  expect  the  behavior  of  the  amplitude  of  a plane  acceleration 

wave  to  be,  at  least  qualitatively,  similar  to  that  of  a purely  longitudinal 

wave.  This  is  for  the  most  part  the  case  but,  as  expected,  in  contrast  to 

the  situation  prevailing  in  the  one -dimensional  case  discussed  in  Ref.l, 

the  coefficients  cd  and  3 , defined  by  (4.34)  and  (4.35),  respectively,  are 
o o 

e 

not  absolute  constants  for  a given  material  and  state  even  though  o,  F and  p 


are  uniform  ahead  of  the  wavefront,  but  vary  with  the  propagation  direction  n. 

Of  course,  once  <S,  F and  pe  are  prescribed  ahead  of  the  wave  then  for  a 

given  n,  r is  determined  by  n through  Eqs.  (3.35).  After  the  unit  vector  r 

has  been  determined  from  Eqs.  (3.35),  the  coefficients  CD  and  3 are  fixed. 

' o o 

Thus,  for  a given  state  ahead  of  the  wavefront,  i.e.,  values  of  F,  <$  and  pe. 


U)  and  3 are  constants  for  a given  n. 

O O 

When  neither  of  the  quantities  ou  ,3  vanishes  then  the  solution  of 

o o 


Eq.  (4.3  3)  is 


CD  t 

a (t)  = \ /l(\  /a  ) - l)e  ° + 1]  , 
o o o ' 


(5.1) 


where 

WPo  ’ (5*2) 

and  a(o)  is  the  value  of  the  amplitude  of  the  wave  at  time  t**0.  It  is 

clear  that  the  behavior  of  the  amplitude  of  a given  plane  wave  is  determined 

by  the  coefficients  cd  and  3 as  well  as  by  the  initial  amplitude  a . In 

o o o 

order  to  discuss  all  possible  cases  which  may  arise  we  first  suppose  that 
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u>  >0  and  3 / 0 , (5.4) 

o o’ 

then,  in  view  of  (5.2),  we  have 

sgn  XQ  - sgn  3q  , (5.4) 

and  then  from  (5.1)  three  distinct  possibilities  arises 

i)  If  sgn  a (0)  » ± sgn  3Q  and  | aQ|  < | Aq| , then  aft)  -»  0 mono- 
tonically  as  t — <*>. 

ii)  If  a (0)  = X(0),  then  a(t)  = a(0). 

iii)  If  sgn  a(0)  = sgn  3q  and  | a (0 ) [ > | \q|  , then  aft)  — <»  mono- 
tonically  within  a finite  time  t^  given  by 

t00  = - (l/u>)&!  [1-  (X  /a(0))]  . (5.5) 

w o o 

We  now  suppose  that 

u>  <0  and  3 i 0 , (5.6) 

o o 7 

then  we  have 

sgn  Ao=-  sgn  3Q  , (5.7) 

and  again  from  (5.1)  three  distinct  possibilities  arise: 

i)*  If  sgn  a(0)=-  sgn  3 and  I X | > I a (0)  I then  a(t)  -‘A 

O ’ o ' o 

monotonically  as  t -»  °°. 

ii)*  If  a(0)  = X then  a(t)  = a(0). 
o’ 

iii)*  If  sgn  a(0)  = sgn  3Q,  then  aft)  — <*>  within  a finite  time  t^ 
given  by 

t.  = - a/Wo)Bn  [1  + I \/aJ  J . (5.8) 

It  is  clear  from  the  foregoing  results  that  the  number  XQ  plays  a funda- 
mental role  in  determining  Whether  the  amplitude  of  an  acceleration  wave 
will  grow  or  decay  as  the  /ave  traverses  the  material.  For  this  reason  we 
follow  the  usual  custom  and  call  XQ  the  critical  amplitude  for  acceleration 

waves  encountering  a homogeneous  steady  state.  We  note  that  if  cu  >0, 

o 7 
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1 

1 


I 

then  the  behavior  of  an  acceleration  wave  propagating  into  a piezoelectric 

semiconductor  which  is  in  a uniform  steady  state  is  precisely  the  same  as 

that  of  an  acceleration  wave  propagating  into  a homogeneously  deformed  ma- 
6 7 

terial  with  memory  * . In  particular,  if  the  initial  amplitude  of  the  wave 
is  less  in  absolute  value  than  the  critical  initial  amplitude,  the  amplitude 
decreases  to  zero  as  the  wave  propagates.  On  the  other  hand,  if  the  initial 
amplitude  is  greater  in  absolute  value  than  the  critical  initial  amplitude, 
the  amplitude  of  the  wave  becomes  unbounded  in  a finite  time.  This,  of 
course,  suggests  the  formation  of  a shock.  As  noted  in  Ref.l,  the  case  ouo < 0, 
which  has  no  mechanical  analogue,  is  the  case  of  primary  interest  and  import- 
ance. The  foregoing  analytical  treatment  shows  that  in  this  case  the  amplitude 
of  the  wave  either  tends  to  Xq  eventually  or  else  becomes  unbounded  in  a 

finite  time.  Furthermore,  note  that  in  this  case  if  a(0)  and  P have  the  same 

’ o 

sign,  the  amplitude  of  the  wave  always  becomes  unbounded  in  a finite  time. 
Moreover, since  a(0)  arises  from  the  thermal  noise,  there  are  always  some  a(0) 
with  the  same  sign  as  (3q. 

Let  us  now  consider  the  behavior  of  a wave  for  which  P vanishes.  In 

o 

our  earlier  treatment  of  one-dimensional  acceleration  waves1  we  noted  that 
the  vanishing  of  Pq  corresponded  to  a linear  material.  Indeed,  while  pQ  also 
vanishes  identically  here  if  the  response  of  the  material  is  linear,  it  may 
also  vanish  because  of  a combination  of  other  factors  even  though  the  response 
of  the  material  is  nonlinear.  For  example,  in  a given  material,  once  <S  and  F 
are  prescribed  it  may  be  possible  to  choose  n [and  hence  r through  Eq.  (3.35)] 
in  such  a way  that  Pq  vanishes.  In  particular  it  is  easily  verified  in  the 
relatively  simple  case  of  a purely  transverse  acceleration  wave  propagating 
in  the  direction  of  the  applied  electric  field  and  of  a principal  axis  of 


homogeneous  deformation  in  an  isotropic  material  that  Pq  vanishes  even  though 
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the  response  of  the  material  is  nonlinear.  In  this  simplest  but  extremely 


important  case  Eq.  (4.33)  has  the  solution 


a (t)  = a e 
o 


-0)  t 
o 


(5.9) 


which  means  that  if  u)^  > 0,  a(t)  is  a monotonically  decreasing  function  of 

time  and  the  amplitude  of  the  wave  decreases  as  the  wave  traverses  the 

material.  On  the  other  hand,  if  cu^  < 0,  a(t)  is  a monotonically  increasing 

function  of  t and  the  amplitude  of  the  wave  increases  without  bound  as  the 

wave  traverses  the  material.  Of  course,  if  in  =0,  then  a(t)  = a so  that 

’ o’  o 

the  wave  propagates  at  constant  amplitude. 

Let  us  recall  from  Eqs. (3.28),  (4.28)  and  (4.34)  that 

e 


U) 


r&A.r.  - v E.r,  ]n.r  . 
2pU  ii  ..  e i i_T  j j 


(5.10) 


In  particular,  Eq. (5. 10)  shows  that  u)q  vanishes  whenever  r is  orthogonal 
to  n.  Thus,  if  vanishes,  either  because  the  response  of  the  material  is 
linear  or  because  <£,  F and  n have  appropriate  values,  then  purely  transverse 
acceleration  waves  will  propagate  at  constant  amplitude  and  not  grow  and, 
of  course,  purely  transverse  shocks  will  not  form.  On  the  other  hand,  in 
this  very  special  case  of  purely  transverse  acceleration  waves,  suppose 
that  3q  does  not  vanish  for  plane  wave  propagation  in  a prescribed  direc- 
tion n,  then  Eq.  (4.33)  reduces  to 


so  that 


!i=Ba2 

6t  “o  » 


a ” a [1-pat] 
o o o 


-1 


(5.11) 


(5.12) 


Of  course,  the  solution  (5,12)  has  the  same  form  as  the  corresponding 

solution  for  acceleration  waves  in  nonheat  conducting  elastic  media.  Note 

that  if  a and  3 have  the  same  sign  then  a shock  will  form  after  a time 
o o 
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t - 1/P  a , 

oo  o O 9 


(5.13) 


and,  as  already  noted,  since  a(0)  arises  from  the  thermal  noise,  a shock 

will  always  form  for  nonzero  3 When  <i>  is  zero. 

o o 

When  vanishes  and  the  acceleration  wave  is  not  purely  transverse, 
it  should  be  clear  from  the  above  discussion  that  the  threshold  condition, 
at  which  the  amplitude  a(t)  just  begins  to  grow,  may  be  defined  by 


U)  - 0 

o 


(5.14) 


where 


o 2pU 


-L  [cA  r ~ V ^ Vi] 


(5.15) 


A most  important  limiting  form  of  (2.16)  [or  (2.7)fe  with  (2.12),  (2.8>7, 


(2.9)3  and  (2.7)5J  is 


(5.16) 


Q 

where  m is  the  mobility  tensor  and  D is  the  diffusivity  tensor,  which  may 
KL  KL 


be  written  in  the  form 


DSo,  = “ > 


(5.17) 


and,  of  course,  we  can  have  m^  = mRL  (E^,  t?^) . In  this  simple  but  important 

limiting  case,  in  which  the  current  is  given  by  Eq.  (5.16)  with  the  mobility 

C 

tensor  m^  and  diffusivity  tensor  constant,  from  (5.15),  (4.18),  (4.11), 
(4.13),  (5.16),  (2.21),  (5.17)  and  the  fact  that  G vanishes  ahead  of  the 
wave  front  we  obtain 


o 2pu 


_i_  v p.eN  JX  .r.  [NTmT  & +u J , 
2PU  ..  e ^ K K.  ll  LLRR  N ’ 


(5.18) 


from  which,  with  (5.14),  we  find  that  the  threshold  relation  is  given  by 


VlrW°- 


(5.19) 


When  the  deformation  is  infinitesimal,  we  have 

F „ 6 „ , £ ■=*  6.E . = E,  U«U=U  , 

rL  rL  ’ R Rj  j R * N o’ 

which  enables  us  to  write  (5.19)  in  the  form 


N m * = Nm  E*  = - U , 
L LR  R L LR  R o’ 


(5.20) 


(5.21) 


\ 


where  it  should  be  recalled  that  N (or  n)  is  the  normal  to  the  plane  wave 
surface.  Equation  (5.21)  is  the  generalization  of  the  well-known  relation 
for  the  threshold  field  obtained  in  Eq. (5.13)  of  Ref.l  for  the  one-dimensional 
case  to  the  arbitrarily  anisotropic  three-dimensional  case  treated  here  for 
the  restricted  limiting  form  (5.16)  of  the  current  equation  (2.16). 


6 . Weak  Waves  and  Induced  Discontinuities 
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Following  Coleman  and  Gurtin  , we  define  a wave  of  order  N as  follows: 


A propagating  singular  surface  £ is  a wave  of  order  N if  the  field  ^ (X, t) 
and  its  first  N-l  partial  derivatives  with  respt  t to  X and  t are  continuous 
everywhere  but  the  Nth  order  partial  derivatives  suffer  jump  discontinuities 
at  £,  but  are  continuous  functions  everywhere  else. 

In  particular,  we  note  the  case  N=2  represents  an  acceleration  wave. 

If  N > 2,  the  wave  is  said  to  be  a weak  wave. 

Our  object  here  is  to  study  the  propagation  and  growth  of  weak  waves  in 
a piezoelectric  semiconducting  material.  It  suffices  to  consider  waves  of 
order  3.  We  confine  our  attention  to  the  study  of  plane  waves  and  we  assume 
that  the  material  ahead  of  the  wavefront  is  in  a state  of  homogeneous  strain, 
is  subject  to  a uniform  electric  field  and  that  the  charge  density  of  the  free 
electronic  fluid  is  uniform  and  constant  prior  to  the  arrival  of  the  wavefront. 

Since  we  are  dealing  with  plane  waves  of  order  3,  from  the  compatibility 


conditions  across  £ we  have 


^rL,Kl-ivr,ucl-brVlt- 

"hVlW 

iKdml\,dm-  Vr*L'  <‘-l> 

lI  ‘ ^r,  KlA  " Wl  < \ ‘ JWtf,  mI  * 

■PU-8*,»»4 

It  follows  from  (2.11)  and  (3.16)  in  essentially  the  same  way  that  (3.19) 

followed  from  (3.18)  that  <5  is  continuous  across  £ and  since  <S  is  con- 

Lj  K L 

tinuous  across  £ also,  from  the  kinematic  condition  of  compatibility  $ is 

L 

continuous  across  E.  From  the  geometric  condition  of  compatibility,  we  have 

iSL,Ki-e\v  A.k’-Kv  *— mi-  <6-2> 

In  a similar  manner  (2.13),  (3.8)  and  (4.15)  imply  that  G is  continuous 

* K,  L 

across  E and  since  G is  continuous  across  E also,  from  the  kinematic  condi- 
tion  of  compatibility  is  continuous  across  E,  but  from  the  geometric 
condition  of  compatibility  and  (2.17),  we  have 

[<Ll  - u5n  , 35-  [N _8J  . (6.3) 

~ Kr^  K ~ R 

Since  for  weak  waves  a^  vanishes, an  immediate  consequence  of  (4.28)  is 


<B)P  - P°V)P»  v ° - ,6-4: 

so  that  the  propagation  condition  for  weak  waves  is  precisely  the  same  as 
that  which  governs  the  propagation  of  acceleration  waves. 

In  order  to  obtain  the  differential  equation  which  governs  the  evolu- 
tionary behavior  of  the  amplitudes  of  weak  waves,  we  differentiate  Eq. (4.1) 
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with  respect  to  t holding  X fixed, and  on  taking  the  jump  across  E in  the 
resulting  equation  and  employing  (3.34),  (4.5),  (6.1)2,  (6.3)  and  recalling 
the  steady  uniform  state  ahead  of  the  wavefront,  we  obtain 

2 - - 2 6l>1 

<WW  - 'oVjp’S  * 0LJ  * *VV  - 2Po°N  ~d’°-  '6-5» 

Similarly,  taking  the  material  time  derivatives  of  (4.10)  and  (4.19),  re- 
spectively, ..aking  the  jumps  across  E in  the  resulting  equations  and  employ- 
ing (2.21)1,  (3.34),  (4.5),  (4.8),  (4.13),  (4.14),  (4.22),  (6.1)  - (6.3)  and 
the  compatibility  conditions  and  recalling  the  steady  uniform  state  ahead 


of  the  wavefront,  we  find  that 


e - cVi  - <V>*\  iVi  • 1 - - Wi  • 

Substituting  from  (6.6)  into  (6.5)  and  employing  (3.36)  and  (4.25),  we 


(6.6) 


obtain 


Let  us  write 


a e n 0 LJ 

(Q.  - P u 6.  )d  -a.  b +2p  u . 

DP  o N ]p  p DP  P Ho  N 6t 


b j " V j > cj  = crj>  l£l=1»  <6-8> 

2 

where  r is  the  unit  eigenvector  of  Q.  corresponding  to  the  eigenvalue  o U . 
~ DP  oN 

If  we  now  contract  (6.7)  with  r^.  and  employ  (3.5),  (3.8)6,  (3.35),  (3.36), 
(4.28)  and  (4.34)  we  find  that  the  amplitude  c of  the  third-order  wave 
satisfies  the  equation 


which  admits  the  solution 


c (t)  - c e 
o 


(6.10) 


It  is  now  clear  that  the  evolutionary  behavior  of  the  amplitude  of  a 
weak  plane  wave  is  somewhat  different  from  that  of  a plane  acceleration  wave 
Which  propagates  in  the  same  direction  at  the  same  speed.  In  particular, 
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the  evolutionary  behavior  of  a weak  wave  is  determined  solely  by  the  sign 
of  u>o  and  is  independent  of  the  initial  value  of  the  amplitude  of  the  wave. 
As  we  have  noted  earlier,  uuq  may  be  negative  in  certain  important  circum- 
stances. When  u>o  is  negative  the  amplitude  of  the  weak  wave  will  increase 
without  bound  and  become  an  acceleration  wave  as  the  wave  traverses  the  ma- 
terial. This  behavior  should  be  contrasted  with  the  manner  in  which  weak 


waves  behave  in  other  media  (see  e.g.,  Ref. 6).  However,  if  the  weak  wave  is 
purely  transverse,  the  amplitude  will  remain  constant  as  the  wave  traverses 


the  material  in  accordance  with  the  relevant  portion  of  the  discussion  in 
Sec. 5. 


Let  us  now  suppose  that  there  exists  a particular  direction  n in  which 
three  real  plane  acceleration  waves  may  propagate.  Let  us  denote  the  ampli- 
tudes, unit  amplitude  vectors  and  speeds  of  propagation  of  these  waves  by 
a<i)  r(i)  and  u * i > i«l  ,2,3,  respectively.  It  follows  from  Eq.  (4.32)  that 
the  amplitude  b ^ of  the  third  order  discontinuity  induced  by  the  accelera- 
tion  wave  of  amplitude  a ^ is  determined  by  the  equations 


(Qjp"PoUN  6jp)bp  dj  ’ 


(6.11) 


d<l,-2p  0a,2r'1>  ♦5i1)r(1>r,1).,1,: 

j j 6t  jp  p jpq  p q 


(6.12) 


where  and  are  given  by  (4.28)  and  (4.29)  with  U replaced  by 

jp  jpq 

0(1>  - B n«l>. 

n n 

Now  suppose  that  the  eigenvectors  of  Q. , (n)  are  distinct  so  that  the 

1 J ~ 

vectors  r^  form  an  orthogonal  triad.  We  may  now  write  b*1*  in  the  form 


bd) . y b(i>r(Qf) 
~ L.  a ~ 


(6.13) 


Since  as  a consequence  of  (4.34)  d(1)  is  orthogonal  to  rf1>,  even  though  D(1) 
is  a root  of  Eq.  (3.37)  Eqs. (6.11)  are  consistent  but  do  not  determine  the 
component  of  b(1)  in  the  direction  of  r(1)  uniquely.  Nevertheless,  when 
Eqs.  (6.11)  are  contracted  successively  with  r^2*  and  r^3*,  respectively, 
and  (3.37)  is  employed,  it  follows  that 

Of  course  Eqs. (6.14)  determine  uniquely  the  components  of  the  induced  dis- 
continuity in  the  two  mutually  orthogonal  directions  Which  are  also  orthogonal 
to  the  direction  of  the  amplitude  vector  of  the  primary  acceleration  wave. 
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